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TITLE 

CRACK  STRESSES:  AN  APPLICATION  OF  IIUSKHELISHVILI  ’  S  EXTHISION 
PRINCIPLE  FOR  MULTIVALUED  MAPPING  FUNCTIONS 


ABSTRACT 

Muskhelishrili's  extension  principle  for  the  case  of  geometries  in¬ 
volving  corners  and  cracks  and  hence  multivalued  mapping  functions  is  con¬ 
sidered.  It  is  shown  that  a  useful  application  of  the  procedure  can  be 
found  in  the  study  of  the  local  stress  distribution  near  the  base  of  a 
crack.  The  structures  of  the  "tress  functions  on  the  Riemann  surface 
defining  the  mapping  function  in  the  neighborhood  of  <rQ  corresponding  to 
the  crack  root  can  be  found  and  the  local  stress  distribution  shown  to 
depend  on  <f>'  (c rQ). 
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I.  INTRODUCTION 


In  the  plane  theory  of  elasticity,  Muskhelishvili ' s*  formulation  of  the 
problem  in  terms  of  a  single  complex  function  defined  in  both  the  occupied 
region  and  its  complement  in  the  plane  is  well-known.  Applications  of  the 
concept  to  regions  defined  by  polynomial  and  rational  mappings  of  several 
basic  regions  have  been  carried  out,  e.g. ,  mappings  of  the  unit  circle  by 
Kartsivadze*.  The  author*  recently  investigated  the  extension  of  the  con¬ 
cept  to  multivalued  mapping  functions  of  the  type  encountered  in  the  de¬ 
scription  of  boundary  corners. 

This  report  is  intended  to  show  a  practical  application  of  the  extension 
principle  when  applied  to  multivalued  mapping  functions.  In  the  study  of 
the  local  stress  distribution  near  the  base  of  a  crack,  it  will  be  shown 
that  the  extension  argument  provides  necessary  information  as  to  the 
character  of  the  stress  function.  A  simple  conversion  of  the  stress  function 
defined  in  the  auxiliary  complex  plane  to  the  stress  concentration  in  the 
physical  plane  is  thus  provided. 


II.  THE  EXTENSION  PRINCIPLE  AND  MULTIVALUED  MAPPINGS 

For  concreteness,  we  consider  regions  in  the  physical  Z-plane  corre¬ 
sponding  to  a  closed  contour  C  and  its  interior  S£.  It  will  be  assumed 
that  a  mapping  relation 


Z  -  coU)  (1) 

exists  which  maps  conformally  the  interior  of  the  unit  circle  in  the  £-plane, 
Sj,  into  s£  and  carries  the  unit  circle  r  into  C.  The  existence  of  corner 
points  in  c  implies  branch  point  type  of  singularities  in  equation  1  lying 
on  the  unit  circle,  e.g.,  consider  the  Schwartz-Christoffel  transformation. 

Muskhelishvili's  original  formulation  depends  on  representing  the 
stresses  and  displacements  in  terms  of  two  complex  functions  <M()  and  \p(i) 
holomorpbic  in  S£  as 


cry  ♦  <rx  -  2RW(t)/a>'(Cn,  (2) 

<ry  -  crx  ♦  2i  rxy  -  2{Ztf£]  [*'(£)/«'(£)] '  +  *'(£}}/«'(£),  (3) 

2jx(u  ♦  iv)  =  I G#£)  -  w(£)<£'(C )/«'(£)  -  ^(£),  U) 

where  primes  denote  differentiation  with  respect  to  £,  bars  denote  complex 
conjugates,  and  n  and  K  are  material  constants.  For  boundary  conditions 
in  terms  of  the  applied  load  on  C,  the  corresponding  conditions  on  the  unit 
circle  £  •  ar  -  eis  can  be  written  as 

4>[ar)  ♦  a>[<r)<f>'  (c r)/u)'  [<r)  *  \p[tr)  ■  g(<r),  (5) 

where  g(c r)  is  prescribed  by  the  mapping  and  the  given  applied  load. 
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The  stress  analysis  requires  the  determination  of  the  two  functions 
4>{  £ )  and  \ //(£)  which  must  be  analytic  in  S£  and  obey  the  condition,  equation 
5,  on  the  unit  circle.  Muskhelishvili rs  extension  principle  consists 
essentially  in  replacing  0(£ )  by  an  extended  definition  of  #£)•  Using 
the  same  form  of  extension  as  Kartzivadze,  we  define  <£(£)  in  the  region 
exterior  to  the  unit  circle,  Sj,  as 

0(C)  =  -cu(C)?'(l/{)/S'(l/C)  -WO,  C«Sj,  (6) 

where 


fu/t)  -  ft  1/C).  (7) 

Clearly  0(£)  as  defined  by  equation  6  is  analytic  in  8£  with  the  possible 
exception  of  the  point  at  infinity.  It  follows  from  equation  6.  that 

0(0  =  -?(l/C)  -  «(1/C>0'(C)A>'(C),  CeS$.  (8) 

It  is  evident  for  the  mapping  class  considered  that  0(C)  as  defined  by 
equation  8  is  analytic  in  S£  with  the  possible  exception  of  C  -  0. 

One  of  the  conditions  on  the  extended  definition  of  0(C)  is  therefore 
the  insurance  of  the  analyticity  of  0(C)  at  C  *  0  in  equation  8.  The 
remaining  condition  arises  from  the  satisfaction  of  the  boundary  condition 
equation  5.  Replacing  0(cr)  in  equation  5  by  the  definition  equation  8,  we 
find 

<f>*Ur)  -  4>~[ar )  ♦  [w*(a)  -  w“(tr)]0,-(cr)/a)'"(cr)  -  g(cr),  (9) 

where  the  notation  f  (cr),  (f+(cr)],  is  defined  as  the  value  of  f ( C )  as 
C-’o- through  Sj,  [Sj],  respectively. 

For  polynomial  mapping  functions,  equation  9  reduces  to  the  so-called 
problem  of  linear  relationship,  i.e. 

0f(cr)  -  0”(ir)  =  g(cr)  (10) 

and  the  solution  can  be  found  in  terms  of  a  polynomial  and  a  Cauchy  integral. 

When  branch  points  in  the  mapping  function  lie  on  the  unit  circle  the 
simplification  equation  10  is  no  longer  possible.  The  mapping  function 
equation  1  must  be  considered  as  defined  on  an  appropriate  Riemann  surface. 
With  little  loss  in  generality,  we  assume  the  appropriate  branch  of  the 
mapping  function  lies  on  the  top  sheet  and  the  branch  cuts  are  intervals  of 
the  unit  circle  joining  the  branch  points.  It  is  evident  that  across  the 
branch  cuts,  a>*[<r)  -  oT(a)  /  0  and  hence  the  simplification  equation  10  is 
no  longer  possible. 

Although  the  effectiveness  of  the  extension  principle  in  the  case  of 
multivalued  mapping  functions  as  a  numerical  procedure  has  little  advantage 
over  other  techniques,  one  can  argue  that  considerable  insight  into  the 
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character  of  the  stress  functions  is  gained.  A  practical  use  of  such  in¬ 
sight  will  now  be  demonstrated  in  the  study  of  the  local  stress  distribu¬ 
tion  near  the  base  of  cracks  occurring  in  geometries  involving  corner 
points. 


III.  AN  APPLICATION  OF  THE  EO'ENSION  PRINCIPLE 

Let  us  assume  that  part  of  C  corresponds  to  a  crack;  thus,  there 
corresponds  an  interval  rc  on  the  unit  circle  which  maps  into  the  crack. 
Let  aQ  denote  the  point  on  rc  which  maps  into  the  tip  of  the  crack  Zq. 
Then 


«*(0  -  ({  -  <r0)Il C),  (11) 

where  1(0  is  analytic  in  S£  and  nonvanishing  in  S£  and  rc.  A  finite 
number  of  branch  points  are  assumed  on  the  unit  circle  necessary  to  the 
description  of  the  corner  points  of  C. 

Since  C  is  a  continuous  closed  curve,  the  branch  cuts  can  be  taken  as 
intervals  of  the  unit  circle.  Furthermore,  it  is  easy  to  show  that  no  cut 
is  necessary  on  rc.  By  a  rotation,  the  crack  can  be  considered  as  falling 
on  the  real  axis.  Then  G>(l/0  is  the  analytic  continuation  of  <w(0  across 
rc  by  Schwartz's  reflection  theorem.  Since  &>(£)  can  be  continued  analyt¬ 
ically  across  rc,  certainly  a>(0  is  analytic  at  £  ■  aQ.  Noting  that 
o>*(cr0)  m  o,  it  follows  that  in  the  neighborhood  of  cr0, 

00 

z  -  Zq  =  Y  -  ^o»n-  (12) 

n  -  2 

If  polar  coordinates  (r, 9)  are  introduced  by 

Z  -  ZQ  =  re1*  (13) 

then  by  reversion  of  series, 

t  -  *0  %  7(Z  -  Zol/So  =  /r/Sg  eW2.  (14) 

From  equations  2  and  3  it  is  clear  that  an  estimate  of  the  stresses 
requires  a  knowledge  of  the  local  behavior  of  the  stress  functions.  For 
simplicity,  we  assume  the  crack  surfaces  are  free  from  applied  load.  Then 
it  is  possible  to  define  g (a)  so  that  g(cr)  -  0  on  rc.  From  equation  9  it 
is  clear  that  0+(cr)  ■  0"(a)  on  rc  for  this  case.  Thus.  d>(£)  as  defined  by 
equation  6.  is  the  analytic  continuation  of  0(0,  USf.  across  the  interval 
(in  the  top  sheet  of  the  Riemann  surface  for  the  mapping  function). In 
particular,  0(0  is  analytic  at  4  •  o0,  thus, 

0(0  •  0'(<ro)  ♦  0'.((ro)(C  -  <r)  + -  (15) 
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is  valid  in  a  neighborhood  of  overlapping  points  in  s£  and  Sj. 
it  follows  easily  that 


ax  *  ay  55  R't 


Zy/Z  [0'(<ro))e“i#/2 


JZm.{o-Q )  /F 


Thus, 


(16) 


To  approximate  equation  3,  the  structure  of  0(£)  is  necessary.  From 
the  preceding  arguments,  0'( £ )/&>'( £ )  is  analytic  at  aQ  except  for  a  simple 
pole,  thus, 


0't£)/a/<£)  =  a.^5  -  c^T1  ♦  aQ  ♦  ax(C  -  aD)  ♦ - .  (17) 

Since  0(5)  and  <*>(£)  are  analytic  at  crQ  then  so  are  $(l/£)  and  &>(l/£);  thus, 
from  equation  B,  </»{£)  has  a  simple  pole  at  aQ.  We  can  write 

0(£)  «  a_j(£  -  ctq)-1  ♦  Aq  ♦  -  crc )  ♦ - .  (18) 

Comparison  of  equations  8  and  18  and  the  series  expansions  above  yields 


A_]_  *  -a_iZ0  -  -0'  (ooJZo/SSg, 


Ao  “  -^(o-o1  "  ao  Zq* 


(19) 


A1  -  0'  (o-0,/°o  ~  s2  a-l/°ro  "  ZqO’ • 


It  is  then  a  straightforward  matter  from  equation  3  to  show 


“  ax  ♦  2i  rxy  % 


✓?e-itf/2 


/rwk.(tr0 ) 


0T^)  0'(^oi  _  4>‘ l*0)*-zi6 

Z 


2<To4 


] 


(20) 


OBSERVATIONS 


It  was  Bhown  in  the  preceding  section  that  the  stresses  in  the 
vicinity  of  an  unloaded  crack  can  be  approximated  by  equations  16.  and  20. 
Since  w'(a0)  can  be  calculated  directly  from  the  mapping  function,  the  key 
computation  is  the  evaluation  of  0'(cro).  Thus,  for  geometries  involving 
corner  points  and  cracks,  the  extension  argument  can  be  used  to  show  that 
the  local  stresses  in  the  crack  vicinity  depend  on  the  evaluation  of 
0'(<?o),  independent  of  the  procedure  adopted,  e.g. ,  polynomial  approxima¬ 
tion,  power  series,  integral  equations,  etc. 
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Generalizations  of  the  arguments  used  in  Section  III  to  a  larger  class 
of  geometries  are  self-evident.  Furthermore,  the  restriction  of  the  loading 
conditions  on  the  crack  surface  can  be  removed  by  a  modification  of  the 
argument.  In  particular,  if  g(o-)  /  0  on  the  interval  corresponding  to  the 
crack  surface^ 


*ltt) 


g(t)dt 

Tt^ry* 


where  (M£)  is  analytic  in  the  neighborhood  of  aD.  Thus,  the  character  of 
0(0  and  hence  0(£)  in  the  vicinity  of  is  known  and  an  argument  similar 
to  that  of  Section  III  can  again  be  carried  out. 
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